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Abstract 

We consider the lattice regularization of A*" = 1 supersymmetric Yang-Mills theory 
with Wilson fermions. This formulation breaks super symmetry at any finite lattice 
spacing; we discuss how Ward identities can be used to define a supersymmetric con- 
tinuum limit, which coincides with the point where the gluino becomes massless. As a 
first step towards the understanding of the zero gluino~mass limit, we present results 
on the quenched low-lying spectrum of SU(2) = 1 Super- Yang-Mills, at /3 = 2.6 on 
aV = 16^ X 32 lattice, in the OZI approximation. Our results, in spite of the quenched 
and OZI approximations, are in remarkable agreement with theoretical predictions in 
the supersymmetric theory, for the states with masses which are not expected to get a 
large contribution from fermion loops. 
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1 Introduction 



Softly broken = 1 supersymmetry gives a natural solution to the long-standing 
hierarchy problem present in the Standard Model, considered as a low-energy effective 
model of a more fundamental theory that includes gravity. In recent years, many 
new results on the non-perturbative regime of super symmetric gauge theories have 
been derived and/or conjectured [Q. However, the only proved results have been 
obtained in the context of extended supersymmetry (N > 1), where the renormalization 
properties are severely constrained Q . From these, rigorous results can also be derived 
in the phenomenologically more relevant cases of = 1 and = 0, through the 
introduction of explicit and perturbative soft-breaking terms Q. Unfortunately, the 
analytical methods fail when the soft-breaking terms become large with respect to the 
scale of the strong dynamics. This is of course the most interesting regime, and the 
lattice is a good non-perturbative method to explore it. 

In this work we will consider the case oi N = 1 Super- Yang-Mills theory (SYM). 
Some preliminary results for gauge observables in this theory have been presented in 
IQ] and confirmed with a different method in The standard lore in SYM was con- 
jectured a long time ago ||6|, ^] and our programme is to confirm this picture on the 
lattice. It is believed that these theories show confinement and a mass gap. Further- 
more SUSY is expected to be unbroken. This expectation is based on the fact that the 
Witten index is not zero [^. As will be explained in this work, SUSY Ward identities 
(Wis) can be used to study the issue of SUSY breaking on the lattice. At low energies, 
the relevant degrees of freedom are expected to be colourless composite fields of gluons 
and gluinos and the lowest-lying spectrum can be described by a low-energy effective 
supersymmetric Lagrangian |^]. The bound states belong to chiral supermultiplets. 
The lowest-lying chiral supermultiplet is expected to be P 

'^ = {n,a,x,F,G), (1) 

where vr = A75A is a pseudoscalar, a = AA is a scalar, and x = ^(^/luF^i/^ is a fermionic 
composite field (here A denotes the gluino field and -F^^y the gluon field tensor). The 
fields F and G (scalar and pseudoscalar respectively) are auxiliary fields with no dy- 
namics. In the supersymmetric limit, these particles should be degenerate in mass Q. 
With supersymmetry softly broken by a gluino mass term, the mass splittings between 
fr, a and x have been computed, to lowest order in the gluino mass, in |Q. 

In this paper we will be concerned with supersymmetric SU{2). Our lattice formu- 
lation, described in section 2, is simply the Wilson action for the bosonic components 
and the Wilson-Dirac operator for the gluinos. The presence of Majorana fermions 
implies a modification of the standard measure of the path integral and the rules of 
Wick contractions. This lattice regularization breaks supersymmetry for two reasons. 
Firstly the breaking of Lorentz symmetry implies that the supersymmetric algebra 
cannot be satisfied at finite lattice spacing. Secondly, the Wilson term induces a finite 
mass for the gluino which splits the supermultiplet. This problem is, however, not 
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very different from the explicit breaking of chiral symmetry in lattice QCD with Wil- 
son fermions and can be handled in the same way The essential point is that the 
symmetry-breaking effects that survive the continuum limit can be eliminated through 
renormalization (if there are no SUSY anomalies), which in this case amounts to a 
proper tuning of the gluino bare mass. This follows from very general principles. Since 
the lattice theory is asymptotically free, the continuum limit can be defined at the UV 
fixed point. In this limit the lattice theory must be equivalent to a theory with the 
same degrees of freedom and all possible interactions of dimension four or less, which 
are compatible with the symmetries of the lattice action. In particular, from the ex- 
act gauge and cubic symmetries of the lattice action, it follows that the only possible 
continuum action is the one corresponding to SYM, augmented by a gluino mass term. 
So even though SUSY is broken at finite lattice spacing in a very complicated way, 
the continuum limit of this theory is softly broken SYM, provided there are no super- 
symmetric anomalies. Through a non-perturbative tuning of the bare gluino mass, a 
supersymmetric limit can be obtained, which coincides with the zero gluino-mass limit 
[10|. In section 3, we review how this can be achieved in practice through the use of 
Ward identities. Moreover, we argue that an important simplification of this method 
consists in the implementation of the OZI approximation. A method to determine the 
order parameter of supersymmetry breaking using WI is also discussed in the same 
section. 

The low-lying spectrum defined in (|l|) can be computed from exponential decay of 
the two-point correlation functions at large time separations, as described in section 4. 
The measurements in the scalar and pseudoscalar channels are completely analogous 
to their QCD counterparts. The fermionic composite x is a new channel that presents 
difficulties, which closely resemble those encountered in glueball spectroscopy in QCD 
and, fortunately, can be overcome with the same type of techniques. 

In sections 5 and 6, we present our spectroscopy results in the quenched and OZI ap- 
proximations. Although SUSY is lost in these approximations, this is a convenient first 
step, which enables lattice calibrations and optimizations of the observables, without 
the high computer cost of dynamical fermion simulations. It is also useful in the under- 
standing of the zero gluino-mass limit in a simpler setting, where the chiral anomaly is 
switched off. As in quenched QCD, PCAC is expected to govern the chiral restoration 
of the theory in the continuum limit. Our results strongly support these expectations. 
Whenever comparison with the results of |p is possible, we find agreement within large 
statistical errors. A preliminary version of our results can be found in |11|. 



2 Lattice gauge theory with Majorana fermions 



We use the lattice regularization of the = 1 SYM, first proposed in ref. |10]. The 
gauge field (the gluon), is the usual link variable Uf^{x) in the fundamental representa- 
tion of the gauge group SU{Nc) (here Nc = 2). The purely gluonic part of the action 
Sg is the usual sum over plaquettes formed by link variables U^{x). The gluino A is in 
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the adjoint representation of the gauge group and satisfies the Major ana condition 

A = = CF, A = F = A^C, (2) 

where the superscript C stands for charge conjugation, the superscript T for transpo- 
sition in spin and colour indices (implicit in the above equations) and C is the charge 
conjugation matrix. The fermionic action is the Wilson-Dirac operator for adjoint 
fermions, i.e. 

x,y x,y 

where M = CA is the fermion matrix 



+ l(mo + 4r)5«''C„^4,,. (4) 

The basic notation and index conventions are given in Appendix As usual, we have 
fixed r = 1. The gauge field V^''{x) is defined by 

V^^''{x) = 2Tt[uI{x)T^U,{x)T% (5) 

with the trace acting on colour indices in the fundamental representation. It satisfies 
the properties |Q] 

= = ]-\ (6) 

where the transposition and inversion act on the adjoint colour indices o, b. Using 
these properties, it can easily be shown that the field V^''{x) transforms in the adjoint 
representation of the colour group. The fermion matrix is antisymmetric: 

M«^(x,y) = -M^^(y,x). (7) 

In the simulations, the fermion fields are rescaled, in standard lattice QCD fashion, by 
a factor ^/{1/2K), with the hopping parameter given hy K = (2mo + 8r)~^. 
The theory is defined by the path integral 

Z[J]= J VUe~^' J PAe-^-f-S--^(^)^(^\ (8) 

where J is an external source of the gluino field A. Note that the fermionic integration 
is only over VX (unlike the QCD case, where we integrate over PAPA). This is a 
consequence of the fact that, with Majorana fermions, the field A is not an independent 
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degree of freedom; c.f. eq. Since the action is at most quadratic in the fermion 
fields, these can be integrated out by the standard change of variables: 

Kix)' = xlix) + lY.4iy)iM~%l{y,x) , (9) 

^ y 

from which we obtain 



where Pf[M] is the Pfaffian of the antisymmetric fermion matrix M (see ref. [Q]). The 
gluino propagator is defined as (using eq. (P)) 

(r{A(x)A(y)}) = {T{\{x)X{y)})C = (11) 
'5^ ln(Z[J])" 



6J{x)6J{y) 



C = {M-Hx,y))C={A-\x,y)), 



where T{...} means time ordering. For the last term of the above equation, we use 
the definition M = CA and the properties of the charge conjugation matrix C; see 
Appendix An important consequence of the Majorana nature of the gluino is that, 
unlike the case of Dirac fermions. Wick contractions of the form A(x)A(y) = M~^(x, y) 
and X{x)X{y) = C^M-^{x,y)C are now allowed. 

Note that our formulation is apparently different from that of [Q], where the action, 
path integral and gluino propagator are obtained in two steps. First, the usual Wilson 
fermion action with Dirac fermions is expressed in terms of Majorana fermions. Then, 
by using a "doubling trick" , which identifies "half the Dirac degrees of freedom with 
Majorana degrees of freedom, all quantities of interest (such as the Majorana fermion 
action, path integral and Green functions) are obtained. With this method (see ref. Q 
for details) the gluino propagator is given by 

{T{X{x)X{y)}) = \[{A-\x,y)) + C-'{{A-\y,x)f)C] (12) 

(here T transposes spin and colour indices). Using the antisymmetry of M and the 
basic properties of C (see Appendix^), it can easily be shown that, on a single gauge 
field configuration, 

A{x,y) = C-^A{y,xfC, (13) 

which implies that the two terms on the r.h.s. of eq. ( |l^ ) are equal. Thus, the gluino 
propagator is identical in the two formalisms (eqs. ([Tl| ) and (|l2|)). Similar identities can 
easily be derived for any correlation function of gluino fields. We therefore conclude 
that the formalism of ref. Q and the one presented here are equivalent; the latter is 
conceptually more straightforward, as it avoids the "doubling trick" . 
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We define the local composite operators that we will use in this work: 



P{x) = A(x)75A(x), 
S{x) = \{x)\{x) 
A^{x) = A(x)7^75A(a;), (14) 

Xaix) = ^P;^(x)(a^^)^^A;^, 
where P^^ is the lattice transcription of the field strength F^j^i, , 

with gQ the bare lattice coupling and Ufj,^ the Clover open plaquette. The reason for 
this choice is that P^i, transforms under parity and time-reversal in the same way as 
F^i, in the continuum; see Appendix Finally, its dual tensor is defined as 

P/iuix) = ^ef,upaPpa{x). (16) 



3 Chiral and supersymmetric Ward identities 



The action described in the previous section is not supersymmetric. However, as 
explained in the introduction, this can be cured in the continuum limit through an 
appropriate tuning of the bare gluino mass (and of operator mixing coefficients, if 
composite operators are involved) [|lO| . This is completely analogous to the problem of 
restoration of non-singlet chiral symmetries in lattice QCD, with Wilson fermions, best 
formulated through Wis |12|. In subsec. 3.1 we review the mechanism of SUSY restora- 
tion on the lattice and show how the gluino mass can be tuned to the supersymmetric 
limit non-perturbatively with the aid of axial and/or SUSY Wis. This procedure may 



be difficult to implement in practice. For this reason, we argue in subsec. |3.2| that the 
tuning of the gluino mass can also be realized in the much simpler framework of axial 



Wis in the OZI approximation. Finally, in subsec. 3.3 we discuss the determination of 



the order parameter of SUSY breaking from SUSY Wis. Unlike the rest of this paper, 
all quantities in this section are in physical (rather than lattice) units. 



3.1 Ward identities and restoration of super symmetry 

For any global symmetry which is broken by the lattice regularization, the bare Wis 
have a form which differs from their continuum counterparts by additional terms arising 
from the explicit breaking. However, the renormalization conditions can be fixed by 
the requirement that the renormalized Wis have the same form as the continuum ones, 
up to terms which vanish in the continuum limit. In our case, there are two global 
symmetries at the classical level: the first is the symmetry under a chiral rotation of 
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the gluino field, which is broken by the anomaly. The second is supersymmetry. On the 
lattice both are broken explicitly. Under axial rotations the following WI is obtained 



V^A^ = 2moP + Xa, (17) 

where P and A^ are defined in eq. (0), niQ is the bare gluino mass and Xa is a 
dimension-five operator with two fermion legs 0. 

All quantities in the above relation are bare and need renormalization. Both A^ 
and P are renormalized multiplicatively, whereas Xa can mix with several operators 
of lower dimension. By using the symmetries of the action, namely parity, gauge 
invariance and cubic invariance, three such operators are identified, in terms of which 
a finite Xa can be defined as 

Xa = Xa + (Za - 1)V^^^ - ZaV^A^ - ZqP^,P^, + 2mP (18) 



Pfj^i, is defined in eq. ( p^ ) and P^,/ in eq. (16). The renormalization condition is chosen 



to be the vanishing of the matrix element of Xa between gluino or physical (colourless) 
states in the continuum limit. Schematically we write lima__>o Xa = 0. This implies 
that the renormalized WI is 

Vf,A^ = 2{mo-m)Zp^P + Q + 0{a), (19) 

where 

Afj, = ZaA^, 
P = ZpP, (20) 
Q = ZqP^^P^^ + ZaV^A^. 

Note that the renormalization of the anomalous term Q is not multiplicative. This fact 
was anticipated in eq. (|l^), where the mixing with V^A^ was artificially split into two 
terms, in order to highlight that one term (proportional to Za — 1) renormalizes the 
current while the other (proportional to Za) renormalizes the operator Pf^uPfiu- This 
WI has the same form as the continuum one, provided we identify the renormalized 
gluino mass with 

rhx = {mo — m)Zp^ . (21) 
For details see the second paper of ref . ||l3[| . 



For the supersymmetric WI, the procedure is analogous |10|. Performing super- 
symmetric transformations on the gluino and gauge fields (see Appendix^, we obtain 
the bare lattice WI 

V^S^ = 2mox + Xs, (22) 

where S^{x) is the lattice version of the supersymmetric current, x is defined in eq. (|14|) 
and Xs is a sum of operators of dimension 11/2. The explicit expressions of and 
Xs has been worked out in Appendix 



^The explicit expression for Xa can be found in [nOj 
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As before, we have to renormalize this relation. The identification of all operators 
of equal or lower dimension, with the same quantum numbers defines a finite (sub- 
tracted) Xs- Imposing the renormalization condition (schematically) lim^^o = 0, 
we obtain the renormalized WI 

V^,Sf, = 2{mo-m)Z-'x + Oia), (23) 

where 



X = Z^x, (24) 
Sfj, = ZsS^ + ZtT^, 

with T^{x) = 7^P"^(x)A"(rE). This has the same form as the continuum WI, provided 
we identify the renormalized gluino mass with 

rhx = {mo - m)Z^^. (25) 

The compatibility of ( p| ) and ( p5| ) implies 

{mo — m)Z^^ = {mo - m)Zp^ . (26) 

Both sides of the equation are functions of mo and, for fixed go, and they vanish at the 
same critical value mo = mcr, as expected. The supersymmetric limit is restored when 

mo = mcr = m{mcr,go) = m{mcr,9o)- (27) 

In principle, the Wis (|l^) and (^) can be used for an explicit non-perturbative 
verification of the above scenario. The supersymmetric WI ( |2^ ) for two independent 
operators Oj, becomes 

^{y,S,{xmO)) + (V^r^(x)O,(0)) = 2^^^^°-^(x(x)O.(0),) (28) 

with i = 1,2. This is valid for x 7^ or for Oj invariant under SUSY transformations 
(otherwise there would also be contact terms). The above two equations (for i = 1,2) 
can be solved for Zs/Zt and (mo — m)/ZT- Similarly, from the axial WI (1£), we 
obtain 

^A_IA(^X/^A,{x)OM) - {P,u{x)P,Ax)OM) = 2 ^"^° -"^) (P(x)O,(0)), (29) 

Zq Zq 

which can be solved for {Za — Za)/Zq and (mo — rh)/ZQ. Thus, at finite UV cutoff 
the ratios of renormalization constants can be determined and, more importantly, 
the tuning of the gluino mass to its critical value can be performed in two independent 
ways. However, since at least two independent amplitudes for each WI insertion are 
needed to tune mo, this might be difficult in practice. 
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Xa Xa 




Figure 1: Examples of the two types of insertion of Xa- (a) OZI diagram; (b) Non-OZI 
diagram. 



3.2 Ward identities in the OZI approximation 

In order to obtain a simple alternative method for the tuning of the gluino mass to 
its critical value, we consider the axial WI in the "OZI approximation". It consists in 
separating the insertions of the operator Xa, appearing in any correlation function, 
into two contributions: those arising from diagrams in which the two fermion fields 
of Xa are not contracted (OZI) and those arising from diagrams in which they are 
(non-OZI); see Fig. ||. Note that the anomaly operator P^yP^u can only arise as a 
counterterm to the non-OZI contributions. Thus, we may render the OZI part of Xa 
finite by defining the subtraction 

Xa,ozi = [Xa + {Z'a — 1)V^j4^ + 2m'P]Q^j (30) 

and imposing lima^o Xa,ozi = 0. A renormalized WI can then be obtained 

V^A'^ = 2{mo-rh'){Z'p)'^P' + 0{a), (31) 

with A'^ = Z'j^Afj^ and P' = Z'pP. This WI is an operator relation, which holds in any 
OZI amplitude. Now, using eq. (|3ll) it is easy to get the point mo = m'(mo,(7o) by 
tuning to zero the ratio 2pozr- 

- (P(x)P(0))|oz. = ^''^ 

However, this is only useful if mo — m' vanishes at the same critical point where 
mo — m = 0, which is not at all obvious. To see how this may come about, we also 
consider the subtractions generated by the non-OZI part of Xa', these are 



Xa 



non—OZI 



Xa- Z'iV^^A^- ZQP^,P^, + 2m"P (33) 

J non—OZI 



The sum of OZI and non-OZI insertions must reproduce the anomalous renormalized 
WI (|T9|). This implies that the OZI tuning of mo — fh' and that of mo — fh occur at 
the same point mcr, provided 

m"{mcr,gl) = Q. (34) 



8 



This is expected to be the case for a theory with more gluino flavours, Nj > 1. In 
such a theory, the previous reasoning holds, but in addition there is an additional 
non-singlet chiral symmetry, which yields the standard non-anomalous PCAC WI. 
The explicit breaking of this symmetry, contained in , can be subtracted through 
renormalization 

X^s ^ + (^/^ _ i)^^Af^s + 2m' P"^', (35) 

where the superscript NS stands for non-singlet current and operator. As usual, the 
renormalization condition is lima^o = 0. Notice that and fh' are identical to 
those appearing in the OZI singlet WI, eq. (^), since the insertions of only have 
OZI contributions. 

In the continuum limit we should recover the expected continuum Wis: 

V,A^' = 2m,P^^ (36) 

V^i^ = 2mxP + Q. 

Notice that, in the continuum, the renormalized gluino mass appearing in both Wis is 
the same. It then follows that 

mo-fh' - fh" mo - m' 
= = (37) 



where = Zp^P^^ . This implies that, at the critical point defined from eq. (^2|), 
eq. ( p^ ) is verified. If we further assume that the dependence of m"{mcr, go) on Nj 
is analytic, this result should also hold in the Nf = 1 case. In this case, the point 
of symmetry restoration is given, up to 0{a) terms, by the vanishing of 2pozi] c.f. 
eq.(ll). 

The computations of the present paper have been performed in the quenched-OZI 
approximation. In this case, the requirement that fh" be analytic in Nf can be relaxed, 
since the OZI two-point correlation functions are identical to the non-singlet ones 
in a theory with A'^^ = 2 (internal fermion loops, which could depend on Nf, are 
absent). Hence, the tuning to the zero gluino-mass limit can be done using eq. (^2|). 
Furthermore, it is believed that, in these approximations, spontaneous chiral symmetry 
breaking occurs because of a non- vanishing gluino condensate [P, |7). If that is the case, 
the OZI pseudoscalar tt should become massless in the zero gluino-mass limit. The two 
requirements m^f and 2p — > give two determinations of mo, which should only 
coincide if chiral symmetry is spontaneously broken. We will see that this is confirmed 
numerically. 



3.3 Order parameter of supersymmetry breaking 

Wis can also be used to study the issue of spontaneous supersymmetry breaking, in 
complete analogy with the study of spontaneous breaking of the non-singlet chiral 
symmetry in QCD described in [|l2|. Consider the supersymmetric WI, 

(V^5^(x)x(0)) = 2(mo - m)Z-\x{x)m) + 5\x){F{x)). (38) 
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The operator F is the renormahzed vacuum expectation value of the bare auxihary 
field of the ^ supermultiplet. Thus it is the order parameter of SUSY breaking. For 
a non-zero gluino mass, we can integrate the previous relation over space-time and 
obtain 

= 2 (mo - m) 5^ Z-\x{x)m) + (^(0))- (39) 

X 

The renormahzed order parameter for spontaneous SUSY breaking can then be com- 
puted as 

(F(0)) = -2 lim (mo-m)^Z-i(x(x)x(0)). (40) 



If this expectation value does not vanish, then SUSY is spontaneously broken. In this 
case the WI also implies that, in the supersymmetric limit, there exists a massless 
fermion field which propagates in the (xx) channel. This is the Goldstone fermion 
of supersymmetry breaking (the goldstino). Hence, another signal of spontaneous 
supersymmetry breaking would be to find a massless fermion in the channel {x{x)x{y)) 
in the limit of zero gluino mass (where this limit is defined by using eq. (p^) or eq. ([29|)). 

4 Correlation functions 

In this section we discuss the method for extracting the lowest-lying masses of the 
pseudoscalar, scalar and spinor particles of the supermultiplet $ (i.e. vf, d and Xi 
respectively), in the OZI approximation. 

The two-point correlation functions used to obtain the tt and a masses are 

Cpp{t) = Y.^P{x,t)P{xoM), 

X 

Cssit) = J2i^{x,t)Sixo,to)). (41) 

X 

In order to compute 2p we need also the following correlation functions: 

CAp{t) = Y.iMx,t)Pixo,to)), 

X 

CpAit) = X](P(x,t)Ao(xo,to)). (42) 

X 

In the computation of all these correlations, the point {xq, to) is held fixed at the origin. 
The correlations Cap and Cpa are antisymmetrized, in order to enhance the signal- 
to- noise ratio (from now on Cap stands for the antisymmetrized correlation function). 
The quantity 2p is obtained from the ratio: 

m - m 

where V^/(x) = [/(x + fl) — f{x — fi)]/2 is a symmetric lattice derivative. 



10 



We now express the two-point correlation functions of eqs. (pll)-(|4^) in terms of 
gluino propagators. The correlation function 

C{t) = Y,{Oi{x)02{xo)) (44) 

X 

of two generic bilocal operators Oi = AFjA, with Fj any Dirac matrix (i = 1, 2), can be 
written as 

C{t)= - ^{Tr[Q{x,t;xo,to;U^)T2Q{xo,to;x,t;Uf,)Ti]) 

X 

+ ^(Tr[Q(f,t;x,t;C/^)ri]Tr[g(fo,to;xo,to;C/^)r2]) (45) 

X 

+ 5](Tr[Q(f, t; xo, to; U^)r2Q{xo, to;x, t; U^)CrJC]), 

X 

where Q{x,t;0,0;U^) = A^^{x,t;(),0;U^) is the fermion propagator, computed on a 
single background field configuration (i.e. the inverse of the fermion matrix). The 
first two terms of the above expression are the familiar one- and two-boundary quark 
diagrams which also appear in the Dirac fermion formulation. The third term arises 
from the extra contractions characteristic of Majorana fermions; see also comments 
after eq. (p!T|). It is also a one-boundary quark diagram. For the cases of interest 
(i.e. eqs. (pl|)-(^3|)), we have that CTJC = — Ti, and the first and third terms are 
equal. In the present work we are interested in the OZI approximation of the above 
expression, for which the second term does not contribute. In other words, we compute 
only one-boundary diagrams. 

For the x mass, we use the two-point correlation function 



x,xo 



= ^/?.(^' t)Qf'p'i^^ ^0, to; u^)pU^o, to) (46) 

This correlation remains unchanged in the OZI approximation, as it only has one 
valence gluino line. For C^{t), we keep to at the origin and vary xq over the whole 
time-slice to- This is obtained in the simulations by solving numerically the equation 

A^";,(f , t; x', t')S^',^(x', t'; xq, to) = 4,^/mo^;.(^> t) (^M-)a/3 • (47) 

The solution of the above is 

Ki3ix,t;xo,to) = Q'^i,{x,t;xo,to;U^)Plil{xo,to){af,^)^,^, (48) 

as can be verified by back substitution. Once T, has been obtained numerically by 
inverting eq. (p7|), the correlation can be easily computed. By using various prop- 
erties (parity, time reversal, etc.) we show in Appendix |^ that the correlation function 
has the form 

[Cx{t)U = Ci{t)6af3 + C2(t)7°/3- (49) 
11 



We recall that the local operator Xa{x) also extends in the time direction, since it 
contains the field strength tensor P^^^. In order to eliminate the contact terms arising 
from adjacent time-slices, we have also considered the operator that has only the spatial 
component of the field strength, defined as follows: 

xf(x) = ^i^^(x)K)„^A^, (50) 
where i,j = 1,2,3. Its two-point correlation function is given by: 

x,xo 

\ i^ij)aa' P^^j{x,t)Q^%>{x,t■,Xo,to)Ptl{xo,to) (^h)^,^ , (51) 

where k,l = 1, 2, 3. 

The fr, a and x masses are extracted from the asymptotic behaviour of the cor- 
responding correlation functions at large time separations. In principle, both form 
factors of eq. ( |4^ ) should decay exponentially in time with the same mass. Hence, we 
can obtain in principle two independent measurements of the x mass. 

5 Particle masses 

We have performed simulations at /? = 2.6. The lattice volume is ^ = 16^ x 32. We 
have generated quenched configurations, separated by 400 thermalizing sweeps, with a 
standard hybrid heat-bath algorithm. We have measured 950 gluino propagators (the 
measurement has been made every 30 sweeps) at four values of the hopping parameter: 
K = 0.174, 0.178, 0.182, 0.184. In H], simulations performed at /3 = 2.3 gave too large 
a value for the a (greater than 1.5 in lattice units). We have opted for a larger /? value, 
hoping to keep all masses smaller than 1 (in lattice units), by reducing the UV cutoff 
a. Statistical errors have been obtained with the jacknife method, by discarding 10 
configurations at a time. Unfortunately the time extension of the lattice seems a little 
too short to isolate the lowest-lying state in the tt channel and two-exponential fits are 
needed in this case. This increases the statistical error in this channel. 

5.1 The TT particle 

We isolate the lowest-lying mass of the correlation function Cpp{t) by fitting it with 
the function 

Cppit) = — exp(-mt) [1 + Za exp(-Amt)l + it ^ T - t) . (52) 
2m 

We have used this two-state fit in order to estimate the lowest mass, by fitting in time 
interval ti < t < T/2 — 1, where T = 32 is the time extension of the lattice. In table |l| 



[ci{t)U = 



12 



K 


ti = Q 


U = 7 


ti = 8 


ti = 9 


U = 10 


0.174 
0.174 


0.869(2) 


0.867(2) 


0.864(5) 
0.876(1) 


0.855(6) 


0.860(7) 


0.178 
0.178 


0.713(2) 


0.711(2) 


0.707(5) 
0.720(1) 


0.698(8) 


0.703(11) 


0.182 
0.182 


0.536(2) 


0.534(3) 


0.531(4) 
0.543(1) 


0.526(11) 


0.528(11) 


0.184 
0.184 


0.433(2) 


0.432(3) 


0.430(4) 
0.441(1) 


0.428(7) 


0.426(9) 



Table 1: The mass of vf as a function of initial fitting time. For each K value, data at the 
top (for 6 < tj < 10) are the results from two-mass fits, whereas the data at the bottom of 
the tj = 8 column are those from one-mass fits. 



we show the results of our two-state fits at the different values of the gluino mass as a 
function of the initial time ti. 

We note that the two-mass fit results decrease with increasing ti in the interval 
6 < ti < 8. This probably signals the weak residual presence of a higher mass state. 
Since, however, from ti = 8 onwards, the values are fluctuating within their error, the 
isolation of two low-lying states is achieved. As ti increases, so does the statistical 
error. We take as the best estimates of the vf mass the results of ti = 8. As a further 
test, we did a one-state fit, also shown in table 2, in the interval 8 < t < 16. The 
results are systematically higher than the ones obtained at the same interval with the 
two-mass fit. This further justifies the necessity for the latter fit. 



5.2 The a particle 

The a particle is noisier than the vr. In table § we show the results of the following 
one-state fit for this correlation: 

Cssit) = Zo + p- [exp (m(t - T)) + exp(-mt)] . (53) 

We take as our best estimates the masses from ti = 10. In all cases, the constant 
Zq was compatible with zero (and orders of magnitude smaller than Zi). This is 
hardly surprising, since, working in the OZI approximation, we are only measuring 
the contribution of the one-boundary diagram to the above correlation function. This 
contribution is equal to the two-point correlation function of the corresponding non- 
singlet scalar operator S^^{x), for which Zq = 0. 

In fig. ^ we present the effective masses for the vr and ct as a function of t, compared 
with the results of eqs. (^) and (^) at K = 0.182. It is clear that the curves fit the data 
well. However, lattices with larger time extensions are needed for the effective mass of 
vf to reach a plateau. In the case of the a, the effective mass plot displays a plateau. 
The straightforward measurement of the a mass in our case must be contrasted to 
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K 


U = l 


= 8 


ti = 9 


ti = 10 


ti = 11 


0.174 


1.13(1) 


1.11(2) 


1.10(2) 


1.09(2) 


1.08(3) 


0.178 


0.99(1) 


0.97(2) 


0.95(2) 


0.95(3) 


0.93(3) 


0.182 


0.83(2) 


0.82(2) 


0.81(3) 


0.81(4) 


0.80(5) 


0.184 


0.75(3) 


0.74(3) 


0.74(5) 


0.76(6) 


0.76(9) 



Table 2: The mass of a as a function of initial fitting time. 



K 


0.174 


0.178 


0.182 


0.184 


2p 


0.379(2) 


0.254(2) 


0.145(1) 


0.096(1) 



Table 3: The values of 2p clS cl function of K. 



the situation in quenched QCD, where the signal for the scalar particle a is extremely 
noisy. 

5.3 The 2p correlation 

The ratio R{t), which gives the estimate of 2/3, is shown in fig. |3| (for K = 0.178). A 
plateau sets in at fairly early times. We perform a weighted average of the results in 
the time interval 9 < t < 14. The 2p estimates thus obtained are shown in table |^ 

5.4 The X particle 

We have obtained the effective mass of x from the spatial correlation function of 
eq. (|5ll). We used all the properties of C^{t) under the discrete symmetries (explained 
in full detail in Appendix B), in order to isolate the two independent components of 
the correlation, namely Ci{t) and C2{t). No effective mass can be obtained naively 
from either Ci{t) or C2{t), owing to large statistical fluctuations. In order to improve 
the signal-to-noise ratio, we implemented the smearing d la APE [^], which is widely 
adopted in glueball spectroscopy. The value of the smearing optimization parameter e 
of ||l^ is fixed at e = 0.5. We have obtained results with two smearing steps, namely 
rig = 4 and rig = 8. In fig. |^ we show the effective mass extracted from Ci{t) and C2{t) 
at K = 0.174 as a function of time (mi and m2 respectively), with = 4 and = 8. 
The results obtained with = 4 show a plateau at small times; this is reminiscent 
of the results for the glueball masses, obtained with the same method. The results 
with Us = 8 have somewhat smaller statistical errors up to times t = 5, but the values 
fluctuate more as a function of t. We plan to study the optimization of the number of 
smearing steps in a future simulation. Since the signal obtained from C2(t) is cleaner 
than that obtained from Ci(t), in this paper we use the former correlation function in 
order to extract the x uiass. We take as our best estimates of the x mass the results 
for 1712 in the interval 3 < t < 5, for = 4, shown in table ^ 
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Figure 2: Effective masses for vf (below) and a (above) as a function of time at = 0.182, 
compared with results of the fits (dotted hues). 



K 


0.174 0.178 0.182 


0.184 




0.93(12) 


0.84(11) 


0.73(10) 


0.66(9) 



Table 4: The values of the x effective mass as a function of K. 



6 The chiral limit 

As explained in section 3, in the quenched OZl approximation, chiral symmetry is not 
broken by the anomaly, but it is expected to be broken spontaneously If this is the 
case, there are two numerically independent determinations of the critical mo^Kc) at 
which the gluino mass vanishes. The first one is obtained upon fitting the vr mass as a 
function of the hopping parameter according to the standard PCAC relation: 

m| = A(l-i-), (54) 

By using this relation we find the chiral limit at 

Kc = 0.18752(9) (55) 
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Figure 3: R{t) as a function of time at = 0.178. The dotted line is its average value in 
the time interval 9-14. 



(with one-mass fits, we obtain Kc = 0.18759(2)). The second possibihty is to extract 
Kc from the relation between two-point correlation functions derived from the WI ( p2[ ) 
in this approximation, which imphes the vanishing of 2p. Extrapolating linearly 2p in 
1/K we obtain 

Kc = 0.18762(4), (56) 

which is compatible with the previous estimate. The behaviour of the OZI component 
of TT is thus similar to that of the vr particle in QCD, and consistent with the expectation 
of spontaneous chiral symmetry breaking. A discrepancy in the two estimates of Kc 
would imply, either that chiral symmetry is not broken (in which case the tt does not 
need to be massless for zero gluino mass, while the WI still has to be satisfied) or, more 
probably, the presence of large 0{a) discretization errors in the correlation functions^. 
We see that none of these possibilities occur. 

In order to obtain the masses of x &iid a at Kc, we have performed a linear ex- 
trapolation in 1/i^ (from the results of the fit for a and from the effective mass for 
x)- The results are shown in fig. ^. It is clear that a and x do not vanish at Kc- This 



*For a recent exhaustive study of discretization effects in QCD, see ref. 14 
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Figure 4: Effective mass of Ci{t) and C2{t) at K = 0.174 as a function of time (mi and m2 
respectively in the text), with = 4 (below) and = 8 (above). 



demonstrates that the model has QCD-hke chiral behaviour, in accordance with the 
expectations of refs. IgI, IToi. The non- vanishing of a at Kc had aheady been seen in 



15|, at /? = 2.3. A comparison with is possible through the ratio of the a mass to 



the lowest glueball mass (0"*""''), which has been measured for both f3 = 2.3 and f3 = 2.6 



in 13] . The resuh is = 1.23(16) at (3 = 2.6 and 1.35(16) at /3 = 2.3. The two 

results are in fair agreement. 

The dependence of the scalar and fermion masses on the gluino mass is compatible 
with linearity. This is expected to be the case in the unquenched theory, where the 
splittings in the multiplet have been computed to lowest order in the gluino mass [^: 

5 

= C + -C'mx, 
6 

= C + lc'mx, (57) 
6 

= C + C'mx. 

C and C' are constants that in general depend on the renormalization scheme PI. We 



^Note that there is disagreement on the results of the splitting between the two papers of ref. H. In 
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Figure 5: Masses of vf, a, x and 2p as a function of 1/K (the mass of vf is squared). The 
extrapolation at Kc is also shown (curves). The value of the 0+"*" glueball mass at /5 = 2.6 is 
also included. 



can compare the ratio of the slopes with the theoretical prediction l^th = 1-17. We 
find 

A = — ^ = 1.3(5 . 58 

dniy. 

Although the error in the slope ratio is quite large, the agreement is encouraging. 
Moreover, the masses of a and x appear to be degenerate at Kc- Of course, this does 
not have to be the case in the quenched-OZI approximation. Nevertheless, the mass 
of the TT particle is expected to change dramatically because of the anomaly, whereas 
there are no reasons to expect a dramatic change of the a and x masses when the 
OZI and quenched approximations are relaxed. In fig. ^ the 0^^ glueball mass at 
(5 = 2.6 is also included. We recall that this particle represents an auxiliary field in 
the $ supermultiplet, and it should be degenerate in mass with the other particles 
belonging to ^ in the super symmetric limit. Its mass appears to be compatible with 
the extrapolated x uiass and 2a below the a mass. This small discrepancy could 
of course be due to the effect of the OZI and/or quenched approximations on the a 



rederiving eq. (|57|), we found agreement with the second paper of 
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mass (in particular, notice that the OZI approximation only affects the a correlation 
function, but neither x the glueball). 

7 Conclusions 

We have measured the low-lying spectrum of SU (2) with one flavour of adjoint mat- 
ter. This is the matter content of supersymmetric SU{2). Our quenched simulation, 
although not yet supersymmetric, is a first step towards the understanding of the fine- 
tuning to the zero gluino-mass limit, which is necessary in the lattice formulation of 
N = 1 supersymmetry with Wilson fermions. We have used the chiral WI to restore 
this limit in this work and discussed how this can also be done in the full theory. 

The numerical results on the spectrum arc quite encouraging. In the full theory 
the theoretical expectation is that SUSY is not broken spontaneously and the low- 
lying states are the components of a composite chiral superfield, i.e. one complex spin 
zero particle and its fcrmionic partner. Of these fields, only the imaginary part of 
the scalar field gets its mass from fcrmion loops, through the chiral anomaly. The 
other states are not expected to change very much when the quenched approximation 
is relaxed. In fact, we have found that these states are approximately degenerate at 
zero gluino mass, already in the quenched approximation. This could be a first signal 
of supersymmetry restoration at zero gluino mass on the lattice. Finally, the vanishing 
of the vf mass at the critical point is an indication of spontaneous chiral symmetry 
breaking, while the finite % mass supports the theoretical expectation that SUSY is 
not spontaneously broken (although this result must be confirmed in the framework of 
unquenched simulations) . 
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A Notation and conventions 

Here we list our conventions. The number of colours is Nc (in this work Nc = 2). The 
number of flavours is Nf (in this work Nf = 1). Latin lower-case indices (of the form 
a,b, . . . = 1, . . . , — 1) stand for colour in the adjoint representation. Latin upper- 
case indices stand for colour in the fundamental representation {A, B, . . . = 1, . . . , Nc). 
Greek lower case indices from the beginning of the alphabet stand for spin (a, /?,... = 
1, . . . , 4) and those from the middle of the alphabet . . = 0, . . . , 3) denote Lorentz 

space-time components (in Euclidean space-time). Latin lower-case indices (of the form 
i,j,k,l = 1,2,3) denote Lorentz spatial components. 
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The gauge field is thus U^^{x), the fermion field is Aq(x), the field-strength tensor 
is F^^{x), etc. It is also customary to write Aq(x) = A^(x)T", F^y{x) = F^j,(x)r", with 
the colour indices in the fundamental representation (e.g. A^'^, etc.) supressed. 

are the group generators, normalized as follows: 

TilT'^T^] = ]^5''\ (59) 

For Nc = 2, these are the standard Pauli matrices T" = a"" /2. 

The Dirac matrices are in the Weyl representation in Euclidean space-time. In 
terms of the three Pauli spin matrices and the unit 2x2 matrix I, they are given by 

We also define 

75 = 71727370, 
1 

C = 7072- 

We often make use of the following property of the charge conjugation matrix: 

C-^ = -C = C^. (62) 
Finally, we take the sign convention eoi23 = 1- 



B Discrete symmetries on the lattice 



In the spirit of ref. ||T^, we list the discrete symmetry transformations of the gauge 
field U^{x) and the fermion fields A(x) and A(x) in Euclidean space-time. The sym- 
metry transformations of the lattice field tensor P^jy(x), the operators xi^)-, x{x) and 
the gauge field in the adjoint representation V^(x) are then easily derived. We also 
derive useful identities between the gluino propagator Q(x; xq; C//i), computed on a 
single background gauge-field configuration, and the one computed on the symmetry- 
transformed background field. The same properties are also derived for the propagator 
of the X operator. Finally, we derive eq. (pS|). 
Under parity V we have: 

(x,to) {x^ ,t^) = {-x,to) 

'ip{x) — > ijj^ {x'^) = ^Q^{x) 

tpix) ^ ip'^ {x'^) = '4){x)-fo 

Uo{x) ^ U^{x'P) = Uo{x) 
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Uk{x) - 


- unx"-) 


= U-k{x) = 


Pok{x) - 


- PS'dx") 


= -P^kix) 


Plk{x) - 


- PKix'') 


= Plk{x) 


X{x) - 


- X'^ix^) 


= lox{x) 


X{x) - 


X^ix''^) 


= X{x)lQ 


Vo{x) - 


- V^ix-^) 


= Vq{x) 


Vk{x) - 




= V_k{x) = 



Ulix - k) 



(63) 



Under time-reversal T we have: 



(x, to) - 


- (x^,t^) 


= -*o) 


■0(x) - 


V'^(x^) 


= Jolbi^ix) 


^p{x) - 


V5^(x^) 


= ^(x)757o 


Uo{x) - 


- Ul{x^) 


= ^7_o(:r) = 


Uk{x) - 


- uUx'') 


= Uk{x) 


Pok{x) - 


~^ Pffkix'^) 


= -Pok{x) 


Pik{x) - 


- PHix^) 


= Plk{x) 


X{x) - 


- X^ix^) 


= 10l5X{x) 


X{x) - 


- X'^ix^) 


= x{x)l-,lo 


^o(x) - 




= Vlo(x) = 


Vk{x) - 


- Vf (x^) 


= Vk{x) 



(64) 



Vo{x-Q) 



These two discrete symmetries result in the following useful properties of the gluino 
propagator Q{x; y; Up) valid on a single gauge-field configuration: 

Q(.x;y;Up) = 7oQ(x^; y^; C/J(2;^))7o, 

Q{x;y;Uf,) = 757oQ(2;^; y^; t^J(a:^))7o75- (65) 

The above properites are identical to the ones valid for Dirac fermions. Also, the usual 
75 property, valid for Dirac fermions, is also true for Majorana fermions: 

Q{x;y;Up) = i->Q{y;x-Up)^-i^. (66) 

We also re-write eq. (p^) in the form 

Q{x-y-Up) = C~^Q{y;x;UpfC. (67) 

We combine eqs. (|66|) and (|67|) to obtain: 

Q{x;y;Up) = C'^l^Q{x-y-UpY^^C. (68) 

The transpose superscript T and the Hermitian transpose one | only act on spin and 
colour indices, i.e. the transposition in space-time is explicit. 
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We now turn to the two-point correlation function of the x field, computed on a 
single background gauge-field configuration, which we denote by [C^(t;U^)]ai3- We 
write 

[C<^{t; Ui,)]af3 = [Qx(^' ^o,to; Uf,)]ap (69) 

x,xo 

in terms of a x-propagator Q^: 

[Qx{x,t;Xo,to;Uf,)]al3 = \ Y1 i'^l^'^)aa' Pliu(.X,t)Q'^^,p,{x,t;Xo,to;Uf,)P^^{xo,to) {(Tpa)f3,/3- 

x,xo 

We shall drop the spinor subscripts a,/3 from now on. All the properties (|65|)-(|6^) are 
also satisfied by Q;^, as can be found by using eqs. (|63[), ( p^ ) and (^). The V and T 
properties of Q^{x,xq; U/^i), listed in eqs. (|6^), imply for the correlation C^{t; U^): 

C^{t; U^) = joC^it; C/J)7o, C^{t; U^) = 757oC^(r - t; C/J)7o75, (71) 

for a single gauge-field configuration with periodic boundary conditions in time. In 
the first of eqs. ( [7l] ) we now average over field configurations and sum over all space. 
This we do by Yl,x,xo I '-'^ l.h.s. and J2-x,-xo I -^^/T r.h.s., in order to 

obtain 

CxW=7oCxW7o. (72) 

Similarly, on the second of eqs. ( [7l| ) we average over field configurations, taking J T>U^ 
on the l.h.s. and / 'DU^ on the r.h.s., in order to obtain 

Cx(t)=757oCx(^-*)7o75. (73) 

The 75 property of eq. (|66[), combined with translational invariance of the propagator 
implies that 

C^{t) = j,C^{T-t)^5, (74) 
whereas from eq. (|67| ) and translational invariance we obtain 

C^{t) = -CC^{T-tfC. (75) 

Finally, combining eq. (|73| ) with (|7^ , and (|7^ ) with (|75|), we obtain 

Cx(0 = 7o4(^)7o, C^{t) = -Cj,C*^{t)j,C. (76) 

Now, let us consider the most general form of the correlation function C^^ as a 
superposition of Dirac matrices: 

c^{t) = j2am, (77) 
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with Ti = /, 75, 7^, 7^75, icT^i/. Using eqs. ([7^ ) and (^), we reduce this expression to 

C^(t) = Ci(t)/ + C2(t)7o. (78) 
Moreover, from eq. ( [73D we derive that 

Ci{t) = Ci(T-t), 

C2{t) = -C2{T-t), (79) 



which, combined with eq. (74), yields that Ci{t) and C2{t) are real. 

By performing analogous manipulations, it is easy to obtain the same results for 
the "spatial" correlation C^{t). 



C Supersymmetric Ward identity 

The supersymmetric lattice transformations for the gluino fields are: 



2 '^/luPfiu^ 



(80) 



where P^y{x), defined in eq. ([15|), has dimension 2. The supersymmetric transforma- 
tions of the link variables are: 



5JJ^,{x) = +iage-f^X{x)U^{x) 
^eUUx) = -iagUUx)e-i^\{x) 



(81) 



In eqs. (^) and (|8l|), e is an infinitesimal fermionic parameter of dimension —1/2. 
These transformations commute with the gauge transformations and give, in the con- 
tinuum limit, the continuum SUSY transformation for the vector supermultiplet. 
The supersymmetric WI is: 



V^S^{x) = moap,P^p,{x)X'{x) + Xs{x). 

The supersymmetric current S^{x) is a dimension 7/2 operator, defined as 

S,{x) = -ap„-fpP''^,{x)V^\x)\\x + /.). 

In the limit o — s- it reduces to the continuum supersymmetric current. 
The operator Xs{x) is given by 



(82) 
(83) 



Xs{x) 



1 



1 



P:^(x)A-(x) - -P^p{x)V^\x)\\x + t,) 



1 



lX-{x-^^)V;i\x-^)P'^p{x) 



(84) 



+ 



1 



-7.75 [\''{x)P;,{x) - P;,(x - t,)V^\x - ^^)\\x) 
g-i^\\x) [a"(x)(7^ - r)faMV;i\x)\'^{x + m) 
23 



where fabc is the structure constant of SU{2). In the naive continuum hmit, Xs is 
proportional to the lattice spacing times a dimension 11/2 composite operator. Xs 
does not contain only terms proportional to r. This is because supersymmetry is 
broken by both the Wilson term and the explicit Lorentz group breaking (due to the 
finiteness of the lattice spacing). The operator Xs{x) is irrelevant at tree level, but 
divergent at one loop, since it mixes with operators of dimension 11/2 or less under 
renormalization ||lO|| . 
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